一种横观各向同性岩体蠕变模型

层状岩体由于层理面或者结构面的存在，在力学上具有横观各向同性的特点，现有的各向同性蠕变模型难以全面反映横观各向同性岩体的蠕变力学特性。为构建能够反映横观各向同性岩体的三维蠕变模型，以能反映瞬时应变、减速蠕变和稳态蠕变特征的黏弹性Burgers模型为基础，基于常泊松比假定，在三维各向同性蠕变本构方程的基础上，按照算子替换的方法，将横观各向同性柔度矩阵代替各向同性柔度矩阵，并考虑了平行和垂直层理方向岩体蠕变力学行为的差异性，推导了横观各向同性岩体的三维蠕变本构方程。根据本构方程的特点，提出了根据平行和垂直方向岩体蠕变试验结果进行三维蠕变本构模型中蠕变参数的辨识方法。将提出的模型应用于三轴蠕变试验参数辨识，从而获得了一套完整的三维蠕变参数，并与试验结果进行对比分析，从而验证了所提模型与方法的合理性和有效性。进一步，指出了传统蠕变试验设计方案的局限性，给出了横观各向同性材料蠕变试验设计建议。研究成果为研究岩体三维蠕变机制提供了新思路，可对岩体蠕变试验设计提供相应的科研支撑。     

南海北部西沙海槽沉积物管状黄铁矿特征及其形成机理*

南海北部西沙海槽S1站位的岩心柱沉积物中广泛发育自生矿物黄铁矿,其形态以管状为主,且具有内部中空的圈层结构。使用扫描电镜、电子探针、LA-ICP-MS、SIMS等测试方法研究了管状黄铁矿的形态及圈层结构,结果显示: (1)管状黄铁矿发育内部中空的圈层结构,其中内圈层(I_py)由莓球状黄铁矿呈五角十二面体紧密堆积组成,外圈层(O_py)由晶形较好晶粒较大的八面体黄铁矿组成,并混有沉积碎屑及钙质生物壳体;(2)内圈层和外圈层分别呈现出贫S富Fe和富S贫Fe的特征,其成因是甲烷渗漏造成的局部还原环境使得As进入黄铁矿中导致晶格空缺或被扭曲,从而促进Ni、Co、Cu、Zn、Pb等微量元素的掺入;(3)内圈层、外圈层发生了明显的硫同位素分馏现象,内圈层中 δ³⁴S 平均为-37.8‰,外圈层中 δ³⁴S 平均为-29.3‰。研究认为,管状黄铁矿作为曾经甲烷渗漏的通道,其生长机制可分为3个阶段: (1)气水通道形成阶段: 向上运移的甲烷流体在沉积物孔隙中逐渐形成气水通道;(2)外圈层形成阶段: 当向上运移的甲烷与硫酸盐发生甲烷厌氧氧化时,逐渐形成晶体较大、晶形较好的八面体黄铁矿外圈层;(3)内圈层形成阶段: 随着甲烷浓度逐渐降低,在气水通道中的微生物作用下,剩余甲烷与向下运移的硫酸盐继续反应形成莓球状黄铁矿内圈层。因此,南海北部的泥岩中大量发育的管状黄铁矿常常与地层中甲烷水合物的存在有关。     

2007-2017年中国沙尘气溶胶的三维分布特征及输送过程

利用最新发布的CALIPSO产品，构建了2007-2017年中国沙尘气溶胶的三维分布，并结合HYSPLIT-4模式和再分析数据，探讨了沙尘的三维输送过程。结果表明：中国的沙尘排放源区主要是塔克拉玛干沙漠和巴丹吉林沙漠，沙尘气溶胶出现频率分别为60%和35%。塔克拉玛干沙漠排放的沙尘主要（50%~70%）停留在源地0~6 000 m高度，少部分向东输送至甘肃和内蒙古；巴丹吉林沙漠排放的沙尘则主要向东输送。中国沙尘排放量在春季最大，向东输送最强；夏季，东亚夏季风限制了沙尘向东输送；秋季，沙尘排放减弱，输送强度和夏季相当；沙尘排放量在冬季最小，输送最弱。夏季，沙尘在输送过程中可被抬升至高度5 000 m以上，春季次之，秋、冬季的沙尘主要在低层大气输送。沙尘在向东输送的过程中被抬升并和当地人为污染物混合变为污染性沙尘，华北地区污染性沙尘出现频率高达30%；输送到海洋的沙尘也会与洋面上（0~3 000 m高度）的海盐气溶胶混合，出现频率约为10%。     

GNSS多天线基线网单历元模糊度同步解算法

提出多天线基线网单历元模糊度同步解算法，其具体实现步骤为：1）根据宽巷组合模糊度易于固定的优点，采用附加已知基线长度约束法同步解算各基线的宽巷模糊度，得到dm~cm级精度的近似基线分量；2）将解算得到的各近似基线分量作为约束，同步解算各基线的基频模糊度，以获取mm级精度的基线矢量。该方法的关键在于检验各历元宽巷模糊度解算的正确性，以获取可靠的近似基线分量，为解算各基线的基频模糊度提供准确的基线先验信息。由于动态情况下各历元观测信息比较少，单纯依赖ratio检验不可靠，提出结合基线误差、单位权中误差、基线网模糊度闭合环及ratio值等对多组宽巷模糊度进行检验，避免ratio值设置不当导致模糊度检验中发生纳伪和弃真问题。实测数据结果证明，该处理方法使得模糊度解算的成功率提高1%~2%，可以获取移动平台更丰富的导航信息，提高其服务能力。     

琼东南盆地华光凹陷天然气水合物稳定带厚度的影响因素

为了探讨琼东南盆地华光凹陷海底天然气水合物稳定带的分布规律，定量研究了静水压力、底水温度、地温梯度和气源组分对水合物稳定带的影响程度。在此基础上，分析了华光凹陷现今甲烷水合物稳定带的厚度分布。最后，综合各因素的历史演化过程，初步探讨了华光凹陷1.05 Ma BP以来天然气水合物稳定带的演化。结果表明:(1)气源组分和海底温度的变化对研究区内水合物稳定带的影响较大；水合物稳定带厚度与海底温度呈良好的线性负相关性。(2)水深超过600 m的海域具备形成天然气水合物的温压条件；超过600 m水深的海域水合物稳定带厚度大部分超过 100 m，其中西北部稳定带的最大厚度超过300 m，是有利的水合物勘探区。(3)华光凹陷1.05 Ma BP以来天然气水合物稳定带厚度经历了快速增厚–窄幅变化–快速减薄和恢复的过程。麻坑群与水合物稳定变化敏感区在空间上具有较好的叠合关系。结合前人的研究成果，推测其形成与天然气水合物的分解释放有关。     

遥感技术在“6·17”丹巴堵江泥石流灾害链灾区应急救援抢险决策中的应用

2020年6月17日,四川丹巴县半扇门镇梅龙沟发生泥石流,阻断小金川河,形成堰塞湖,造成重大财产损失与人员伤亡。为掌握第一手灾情,辅助开展应急救援,以国产高分二号卫星数据为主要数据源,结合无人机航空摄影数据和灾后现场调查等资料,开展了“6·17”丹巴堵江泥石流灾害链灾区应急调查与分析。研究认为: 灾区影响人口约6.5万人; 国网丹巴县供电公司应急救援队与小金县联通通信应急抢险救援队距离最近,半扇门中学等5所学校可作为临时安置避难场所备选; 受影响矿山企业6家、重要水库1座、水电站1座; 泥石流沟总体呈NW—SE向展布,掩没面积约14.27万m²; 堰塞湖湖面面积约1.03 km²,淹没区约49.8万m²,路基路面全毁路段约2.1 km,下游疑似灾害隐患点9处,因灾受损民居25处、桥梁5处; 优选了2条灾后通达性较好的通往泥石流重灾区的救援生命线。借助遥感技术深入分析泥石流灾害链影响人口等灾情先期研判、灾情灾损应急调查及应急救援生命通道优选,对我国西南山区类似条件地质灾害隐患点的应急抢险救援、防灾减灾工作具有重要的指导意义。     

湍流度随高度的变化及其对城市宏观地形的依赖

分别构建广州主建成区垂直比例尺为1﹕2 000、1﹕1 000和1﹕500的3个建筑物模型，利用大型边界层风洞，在西北和东南两风向下，基于中性流模拟分析了复杂城市地形下湍流度随高度的变化及其对宏观地形的依赖。结果表明：风廓线指数α与不同高度的湍流度之间的关系密切，利用现有模型，根据4类粗糙度边界层和不同垂直比例尺，可确定相应的湍流度随高度变化模型的主要系数，预测精度高。城市地形下最大湍流度面发育在0~0.2 h之间狭窄的范围内。用湍流度形态指数β来表征湍流度随高度的变化，无论城市屋脊还是平坦地形，随着风程区的延伸，廓线的指数α升高，湍流度形态指数β降低。表明同一高度湍流度值具有由迎风区、丘顶区向背风区增高，沿风程逐渐增大的规律，对地形部位和风程的依赖性强，与来流翻越简单地形时的特征一致。     

Progress in Danube Loess and Paleoenvironment Study

Danube loess, located at the westernmost part of the widespread Eurasian loess belt, is one of the most well-preserved aeolian dust deposits in Europe with the basal age dating back over 1 Ma. Owing to its unique location and formation processes, Danube loess plays an important role in understanding the linkage of climate changes in different climate regimes over Eurasian continent and paleoclimatic changes in middle-high northern latitudes. Major research advances of the Danube loess, including compositions, provenance, chronology, stratigraphy and paleoenvironment changes were systematically reviewed and the focuses of future research were suggested. To better understand loess provenance, the compositions of sediments in all the potential source areas should be investigated, and methodological study of provenance indicators and application of multi-proxies approaches need to be carried out. Chronologically and stratigraphically, feldspar luminescence dating and tephrochronology methods should be widely used, and the establishment of a uniform stratigraphic framework and correlation scheme should be constrained by precise chronology. In terms of the paleoenvironment, it is necessary to conduct high-resolution multi-proxies reconstruction of paleoenviroment for the entire loess-paleosol sequence, extending from the last interglacial period (130 ka to present) to interval of last 1 Ma. The future progress in Danube loess is crucial to improving the study of continental paleoclimate comparison and paleoclimate dynamics over the Eurasian loess region.     

西藏甲玛铜多金属矿床磁黄铁矿标型矿物学特征及其地质意义

甲玛矿床位于冈底斯成矿带东段,是西藏地区最大的铜多金属矿床之一。磁黄铁矿是甲玛矿床最常见的金属矿物之一,其标型特征不仅反映其自身形成环境,对其形成机制和矿床成因也具有指示意义。文章选取产于不同岩性中的磁黄铁矿矿石样品,利用矿相学、X射线衍射和电子探针分析等手段对磁黄铁矿的形态、成分和结构进行了分析研究。研究表明,甲玛矿床的磁黄铁矿主要分布在距离岩体中心较远的矿区远端矽卡岩和角岩中。磁黄铁矿的晶胞参数和粉晶X射线衍射曲线显示矽卡岩中的磁黄铁矿主要为高温六方磁黄铁矿,角岩中的磁黄铁矿为高温六方磁黄铁矿和低温单斜磁黄铁矿的交生体,但主要以低温单斜磁黄铁矿为主。通过对矽卡岩和角岩中的磁黄铁矿进行电子探针测试,结果显示:矽卡岩中的磁黄铁矿中w(Fe)为60.09%～60.71%,平均为60.38%,w(S)为38.18%～38.69%,平均38.35%,化学分子式为Fe_8S_9～Fe_(10)S_(11);角岩中的磁黄铁矿中w(Fe)为59.05%～59.57%,平均为59.10%,w(S)为39.28%～39.95%,平均39.59%,化学分子式为Fe_5S_6～Fe_7S_8。根据以上矿物学特征,笔者进一步探讨了该矿床磁黄铁矿的沉淀机制:炽热的岩浆热液上涌,与碳酸盐岩地层和碎屑岩地层接触发生相互作用,并有大气水的加入,使得成矿流体在角岩中先快速降温,形成高温六方磁黄铁矿和低温单斜磁黄铁矿的交生体。同时,大量的含矿热液形成,并充填于有利的成矿空间(主要为层间破碎带)沉淀成矿,形成矽卡岩矿体,然后流体在矽卡岩矿段中经历缓慢降温,形成高温六方磁黄铁矿。结合矿床地质特征和相关元素地球化学特征,认为甲玛矿床类型为斑岩-矽卡岩型。     

利用综合评分方法评估CMIP5模式对北极海冰的模拟能力

北极海冰冰盖自20世纪以来经历了前所未有的缩减，这使得北极海冰异常对大气环流的反馈作用日益显现。尽管目前的气候模式模拟北极海冰均为减少的趋势，但各模式间仍然存在较大的分散性。为了评估模式对于北极海冰变化及其气候效应的模拟能力，我们将海冰线性趋势和年际异常两者结合起来构造了一种合理的衡量指标。我们还强调巴伦支与卡拉海的重要性，因为前人研究证明此区域海冰异常是近年来影响大尺度大气环流变异的关键因子。根据我们设定的标准，CMIP5模式对海冰的模拟可被归为三种类型。这三组多模式集合平均之间存在巨大的差异，验证了这种分组方法的合理性。此外，我们还进一步探讨了造成模式海冰模拟能力差别的潜在物理因子。结果表明模式所采用的臭氧资料集对海冰模拟能力有显著的影响。